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Abstract. Kapitza-Dirac-Talbot-Lau interferometry (KDTLI) has recently been 
established for demonstrating the quantum wave nature of large molecules. A phase 
space treatment permits us to derive closed equations for the near-field interference 
pattern, as well as for the Moire-type pattern that would arise if the molecules 
were to be treated as classical particles. The model provides a simple and elegant 
way to account for the molecular phase shifts related to the optical dipole potential 
as well as for the incoherent effect of photon absorption at the second grating. 
We present experimental results for different molecular masses, polarizabilities and 
absorption cross sections using fullerenes and fluorofullerenes and discuss the alignment 
requirements. Our results with Ceo and C70, C60F36 and C60F48 verify the theoretical 
description to a high degree of precision. 
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1. Introduction 

The quantum wave nature of matter has become a corner stone of physics over many 
decades, and current interest in de Broglie interferometry with electrons [D [21 G], 
neutrons [I], atoms [SHE], and molecules [HE] ranges from demonstrating fundamental 
quantum phenomena to advanced applications in the materials sciences and in quantum 
metrology. All these experiments require optical elements for the coherent manipulation 
of matter waves. While clean solid surfaces and bulk crystal structures are well-adapted 
to the diffraction of electrons and neutrons with de Broglie wavelengths in the range of 
1.. 1000 pm, it is often necessary to tailor the beam splitters, lenses, and wave guides to 
the specific particle properties in atomic and molecular applications. 

For complex molecules, nanofabricated gratings were demonstrated to act as beam 
splitters for far- field diffraction [9] and near-field interferometry [TO] . However, these 
experiments pointed already to the importance of van der Waals interactions between 
the molecules and the diffraction grating, which largely exceeds the effect observed with 
atoms [UJ because of the high molecular polarizability and their comparatively low 
velocity. The interaction time with a 500 nm thick grating amounts to only 5 ns at a 
beam velocity of 100 m/s, and yet the matter wave phase shift can attain the value 
of several radians in the center of the slit opening. The interaction effect gets even 
stronger close to the slit walls, to a degree that the wave front distortion can no longer 
be described by a phase shift alone [12J. For particles with increasing polarizability 
this strong influence of the grating interaction leads to prohibitive requirements on the 
velocity, as discussed in [13]. It is therefore appealing to replace material gratings by 
structures made of light, which offer the additional advantage of being indestructible, 
highly transparent, and easy to tune and modulate. 

Bragg diffraction of free electrons at a standing light field was already proposed by 
Kapitza and Dirac in 1933 [Hj, but nearly seventy years passed before the idea was 
experimentally implemented [T5]. In contrast to that, the first optical phase grating 
for atoms was already realized in 1983 [T6l [TT] when a standing laser light field was 
tuned near to an atomic resonance in order to perform Raman-Nath ('thin grating') 
diffraction of a supersonic sodium beam. A related investigation [TS] then focused on 
atomic diffraction in the Bragg regime ('thick grating'). These ideas were later extended 
to atom diffraction [19J and interferometry [201 EI], also to the time domain [221 [231 [21] , 
and to the manipulation of Bose-Einstein condensates [251 EE] • 

The working principle of all phase grating examples is the same: A coherent laser 
beam creates a periodic pattern of the electrical field. This couples to the particle's 
polarizability, shifts the energy and thus imprints a phase pattern on the transmitted 
matter wave beam. Its evolution into a modulated particle density distribution can then 
be observed further downstream. 

Large, hot molecules in thermal beams often exhibit broad absorption lines. Light 
will therefore mainly couple in a non-resonant fashion. However, for most molecules one 
can still find a suitable range of wavelengths where the light-molecule coupling allows 
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Figure 1. Setup of the Kapitza-Dirac-Talbot-Lau interferometer: it consists of two 
material nanostructures (Gl, G3) and a standing light wave (G2). The latter is 
realized by a cylindrical lens focusing a 532 nm laser beam onto a mirror. The three 
structures (G1,G2,G3) have the same period of 266 nm and are separated equidistantly 
by 105 mm. For detection, the third grating (G3) is shifted in small steps over the 
molecular interference pattern. The transmitted molecules are detected in a quadrupolc 
mass spectrometer. 

one to imprint a local matter wave phase shift of the order of A$ = n. 

The first application of optical phase gratings to large, hot molecules was 
demonstrated with Ceo in a far-field diffraction experiment [27j. The combination with 
near-field diffraction was suggested in [28J and recently implemented in a Kapitza-Dirac- 
Talbot-Lau interferometer (KDTLI) [13] , as shown in Fig. [TJ 

The general idea behind the KDTLI design has been described elsewhere [2"8"| [TB"] : 
A first absorptive mechanical structure, Gl, (in the present experiment d = 266 nm) 
prepares the required spatial coherence for illuminating the optical phase grating, G2. 
Quantum interference then explains the appearance of an approximate self-image of 
G2 at the position of the third mask G3. This molecular density pattern is scanned by 
shifting G3 across the beam while counting the transmitted molecules as a function of the 
mask position. This scheme was exploited to perform quantum interference experiments 
up to extended polyatomic molecular chains [13] and to determine electrical molecular 
properties [2Uf I3T)]. 

In the present article we now provide a refined theoretical description of Kapitza- 
Dirac-Talbot-Lau interference, putting special emphasis on the proper incorporation of 
the influence of photon absorption in the second grating. We give a closed expression 
for both the quantum interference visibility and the fringe contrast one would expect 
if a classical Moire description were correct. We then compare this to the measured 
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interference curves of Ceo and C70 which are in nearly perfect agreement with the 
quantum result. We also apply the KDTL concept to studying the fluorofullerenes 
C60F36 and C60F48. This allows us to assess the influence of mass, absorption cross 
section, and optical polarizability on the interference of large particles. In comparison to 
earlier Talbot-Lau experiments [31] . the new Kapitza-Dirac-Talbot-Lau interferometer 
now also allows us to establish a significantly improved fringe contrast. 

2. Theory of the Kapitza-Dirac-Talbot-Lau interferometer 

The Kapitza-Dirac-Talbot-Lau interferometer is a derivative of the standard Talbot-Lau 
interferometer [32J [33], El ID] , obtained by substituting the central grating mask with the 
optical phase grating created by a standing laser beam. In the simplest configuration 
all gratings, material and optical, have the same grating period d, given by one half of 
the laser wavelength, d = Al/2. The passage of the matter wave beam through each 
grating may thus transfer integer multiples of the grating momentum p^, = h/d = 2h/\i 
to the transverse motion in the beam. These different diffraction orders interfere further 
downstream, leading to a resonant enhancement at integer multiples of the Talbot length 
Lt = cP/AdB, which is determined by the de Broglie wave length AdB of the molecules 
(which ranges between 1pm and 5pm in our experiment). The emerging interference 
pattern at the position of the third grating thus displays a strong dependence of the 
interference fringe visibility on AdB, as determined by the longitudinal velocity v z of the 
beam. 

Unlike with material gratings, in the KDTLI we must also consider the possibility 
that one or more laser photons are scattered or absorbed while the molecule traverses 
the standing light wave. The associated incoherent transfer of transverse momentum 
may strongly blur the fringe pattern. If the absorption is followed by an immediate 
isotropic reemission process, the transverse momentum shift may take any value up to 
the photon momentum H/\l- However, in many large molecules the absorbed photon 
energy gets stored for a rather long time, either in metastable excited states or, after 
rapid internal conversion, in the vibrational degrees of freedom, which do not decay 
over the time scale of the experiment. The associated transverse momentum transfer is 
then an integer multiple of the photon momentum h/X^, corresponding to one half of 
the grating momentum pd- An odd number of net photon momenta will thus kick the 
molecular wave such that the fringes get blurred maximally, while an even number will 
have a much weaker effect. We note that related physics has already been described in 
the context of far-field diffraction [33] and Mach-Zehnder interferometry [35J with atoms 
before. 

In order to describe the interplay of coherent diffraction and the incoherent 
effect of photon absorption we follow the phase space formulation of Talbot-Lau 
interference based on the Wigner function, as presented in [36]. It provides a transparent 
representation of all relevant coherent and incoherent phenomena, and it permits us to 
calculate the quantum interference pattern on an equal footing with the possible moire- 
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type structures, which might arise already due to classical mechanics in this setup. This 
comparison with the classical description is required if one wants to establish that the 
observed fringe pattern in a molecule interference experiment is caused by a genuine 
quantum interference effect. 



2.1. The light-grating interaction 

We start by collecting the necessary ingredients for describing the effects of a light 
grating on the motion of a beam of polarizable particles. Taking the direction of the 
particle beam as the z-axis, we set the retro-reflected, basic Gaussian laser mode in the 
perpendicular x-direction. The time averaged intensity of the standing light wave is 
then given by 

I\?,y,z)= exp sin vr- , (1) 

7TWyW z \ Wy w z z ) \ a J 

where w y and w z denote the laser beam waists in the vertical and the longitudinal 
direction, and P is the laser power. In the following, we assume the particle beam 
height to be sufficiently small compared to w y , such that the dependence on the vertical 
y-direction can be safely neglected. This is approximately the case in our experiment, 
see below. 

The standing light field will in general induce dispersive and absorptive forces on 
a molecule. The first type, due to the conservative optical dipole force, is described by 
the potential 

V{x,z) = -^l(x,0,z), (2) 

where is the real part of the polarizability of the particle at the laser frequency 
uj = 2ttc/\l (related to the polarizability in SI units by a§i (uj) = AnEQaS). 

Treating the effect of the grating potential in the eikonal approximation, a 
traversing quantum wave acquires a position dependent phase shift which is calculated 
by integrating the potential along a straight line, 



1 f°° ( x\ 

00) = V(x,v z t)dt = o sin 2 (jr-J. (3) 

The maximal shift 

00 = 8^^ (4) 

hC WyV Z 

thus increases linearly with the optical polarizability and laser power, and it is inversely 
proportional to the molecule velocity v z . The justification of this eikonal approximation 
and its range of applicability are discussed in some detail in [12]. As shown there, it 
is well justified for the molecular masses and polarizabilities accessible with the current 
experimental setup. 

The second type of momentum exchange between the light and the molecules is 
the radiation pressure force due to photon absorption. The molecules used in the 
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present experiment are sufficiently large and internally complex that it is justified to 
both ignore any reemission and to take the absorption cross section constant even after 
the absorption of several photons. In this case all absorption events can be described as 
being independent and as only determined by the absorption cross section cr a b s at the 
laser frequency to = lisv. This cross section can often be related to the imaginary part 
of the polarizability using Mie theory, cr a b s = 4:iiu/c x Im (a (a;)), but we will treat cr a b s 
as an independent parameter in the following. 

The photon absorption rate is determined by the incident photon flux I(x,z)/hu 
and a abs , 

V(x,z) = j^I(x,z) . (5) 

Below, the radiation pressure effect on the molecular beam will be described by the 
position dependent mean number of absorbed photons. The latter is obtained from 
the photon absorption rate, in analogy to the eikonal approximation, by a straight 
integration along the longitudinal motion of the molecule, 



n(x) = I r(x,v z t)dt = no sin 2 (71-— j. (6) 

The maximum mean number of absorbed photons n is found in the anti-nodes of the 
standing light wave, and it is given by 

8 0"absAL P . . 

n = —~= — : . (7) 

V27T he w y v z 

The values of O an d n defined in (J3J) and (J7|) are the two key parameters describing the 
molecule-light interaction, and they will appear in the closed formula for the quantum 
interference visibility below. 

Classical description. We note that the momentum- changing effect of photon 
absorptions does not differ in the quantum and the classical description of the molecular 
motion. On the other hand, the classical effect of the optical dipole force due to the 
potential ^ should be treated in analogy to the eikonal approximation of the quantum 
case. This is done in terms of the momentum kick Q(x) obtained by integrating the 
dipole force along the same straight line as in the eikonal treatment [36J, 

^ / x f°° dV , , , irh , / x\ 

Q(x) = - J (x, v z t) dt = —0o sin [2ir-J , (8) 

with 0o from (jlj) . (The h in the prefactor cancels Planck's constant in O rendering the 
equation classical.) 

As we will see below, the spatially periodic focusing of classical particles due to the 
dipole potential (jHJ) may result in a regular molecule pattern behind the light grating, 
though distinctly different from the quantum prediction. To perform the quantum and 
the classical calculations it is useful to formulate the effect of the grating passage in a 
common framework, the Wigner-Weyl phase space representation. 
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2.2. Phase space formulation of the light-grating transformation 

The most important part in describing the KDTL-interference is the transformation of 
the particle beam state as it passes the second grating. We consider the Wigner function 

"frlO = A / dse^(x-^\p\x + ^ (9) 



2ixh J " 2 iri 2' 

of the transverse quantum state of motion p of the molecular beam [371 [38j [36], where 
x and p denote the position and momentum coordinates in phase space, and we first 
assume the longitudinal velocity of the molecules to be given by a definite value v z . 

After passing an arbitrary grating in eikonal approximation, the transformed beam 
state can always be written as [36] 

w'(x,p) = dp T(x,p-p )w(x,p ) . (10) 



That is, the momentum dependence of the quasi phase space distribution gets modified 
by a convolution, while its position dependence is at most affected by a multiplication. 

Let us now discuss the grating transformation for an arbitrary, <i-per iodic 
distribution of the light intensity. We first assume that there is no absorption, 
Cabs = 0, so that the grating transformation is entirely coherent. The phase shift 
(f) (x) then relates the wave function ip in front of the grating to the one behind, 
ip' (x) = exp (i(p (x)) ip (%)■ In phase space representation this coherent transformation 
is described by the convolution kernel [3S] 

T coh (x,p) = ^J dse^exp [z0 (x - |) - *0 (x + £)] . (11) 

Noting the periodicity of <p(x) we define the Fourier coefficients 
i rd/2 

bj = - exp(t<f)(x))e- 2nijx/d dx } (12) 

« J-d/2 

so that the coherent kernel takes the form 

T coh (X,p) = Yl h i h ]-m eX P ( 2nim ~^) 6 (P~ - Pdj ■ (13) 

This is a periodic comb of delta-functions separated by integer multiples of the grating 
momentum p^ = h/d. It serves to populate the different diffraction orders in (1101) as 
the quantum wave passes the grating. 

In a second limiting case, we now consider the grating transformation for a vanishing 
dipole force, = 0, but maintain a finite absorption cross section. It is reasonable to 
assume that the final detection efficiency of the beam particles is practically independent 
of the number of absorbed photons. In this case, their motional state gets effectively 
replaced by a statistical mixture whose components differ by momentum translations of 
integer multiples of the photon momentum H/Xl — Pd/2- These multiples correspond 
to the difference in the number of photons absorbed from the left and from the right 
side. We denote as Prob (k; x) the position-dependent probability distribution for the 
exchange of k e Z net photon momenta. The mixture can then be written in phase space 
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representation as w'(x,p) = J2k P r °b (k; x) w (x,p — kpd/2). That is, the statistical 
redistribution of the momenta due to photon absorption is described by the incoherent 
kernel 

T ahs (x,p) = Y,^(k;x)5(p-k^y (14) 

Expanding the periodic position dependence of Prob (k; x) in a Fourier series, 

Prob (jfe; x) = J2 P f ) e 2nijx/d J (15) 

it takes the form 

T abs (x,p) = ^ xld Pfh (p - k?f) . (16) 

Let us now specify the probability distribution (I15p in terms of the mean number 
of absorbed photons n(x). Since the absorption events are taken to be statistically 
independent the probability Prob (k; x) for the net gain of k photon momenta is given 
by 

00 -n( \ 

Prob (Jfe; x) = Prob (%) e ~ n(X) - (17) 

n=0 n ' 

Here, n (x) is the position-dependent mean number of photon absorptions characterizing 
the Poissonian distribution and Prob (k\n) is the probability for the net transfer of k 
photon momenta towards one side, conditioned on the absorption of exactly n photons. 
Since absorptions from the left and from the right occur with the same probability in a 
standing light wave, the latter is given by the distribution of a one-dimensional, balanced 
random walk with n steps, 

1 I (^."wo) if n + k even , , 

Prob(fc|n)= — <^ Ufc+n)/ 2 ; ^g) 

otherwise. 



In order to perform the average over the Poissonian photon distribution in (IT7I) we first 
calculate the characteristic function of (|18p by means of the binomial theorem, 

X(£|n) = ^exp(-27ri^)Prob(A;|n) = [cos (2^)]" . (19) 

The characteristic function of the averaged distribution fTTTj) thus takes the simple form 

_ _ / . » v^v \n (x) cos (27r01 n ... 

X ;i = Y,-^- 2 — r exp -n x 
A — ' n! 

n=0 

= exp{-n(x) [1 - cos(2tt0]} • (20) 

The inverse Fourier transform of ( 1201) yields the required probability fTTTj) in terms of 
the mean number of absorbed photons, 

Prob (k; x) = exp (— n (x)) Ik (n (x)) , (21) 

where the /„ (x) are modified Bessel functions of the first kind. 
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So far, the integral kernels T co h and T a b s , which describe the coherent and the 
incoherent part of the light grating interaction, were discussed separately, see Eq. ([TBI 
and ( ITB"]) . For realistic molecules the dispersive and the dissipative light forces coexist, 
and their combined contribution is described in the eikonal approximation by a single 
transformation (flOl) . whose kernel is given by the convolution of T co h and T a b s , 

T (x, p) = J dq T coh (x, p-q) T abs (x, q) . (22) 

Classical description. As an advantage of the phase space formulation, it is easy to 
describe in the same framework how the molecules would move if they were classical 
particles. One merely replaces the Wigner function by the classical phase space 
distribution, which is a proper probability density. Both the quantum Wigner function 
and the classical phase space distribution experience the same shearing transformation 
as they evolve freely between the optical elements, see the discussion in Also the 
passage through a grating can be expressed in the same form ( TTDT) in both cases, though 
the integral kernels differ of course. The classical kernel due to the dipole force takes 
the form (P2J) 

T cl (x,p) =5(p-Q(x)), (23) 

where Q(x) is the classical momentum kick of Eq. (jSj). The effect of a photon absorption, 
on the other hand, is described by the same kernel (1T41) as in the quantum case, since 
it effects the same momentum change on the motional state, irrespective whether the 
center-of-mass motion is described by classical or quantum dynamics. To obtain the 
combined effect of dispersive and absorptive light forces, one can again concatenate the 
two transformations. 



2.3. Evaluating the Kapitza-Dirac- Talbot-Lau effect 

We are now in the position to calculate the interference pattern expected for the KDTLI 
in the same spirit as it was done for purely coherent grating interactions in [36]. The 
final beam state is obtained by applying the appropriate sequence of free evolution and 
grating transformations to its Wigner function. Starting with a spatially completely 
incoherent but monochromatic beam in front of the first grating one thus obtains the 
spatial density distribution in front of the third grating by a final integration over the 
momentum variable, 



2J / dx dpTi(x )T coh ix L, 2p — p z - k— J 

lc» J \ Pz J 



W 3 [X) (X 

x Prob Hz, x — ■ ( 24 ) 

Here L is the distance between the gratings and Ti(xq) G {0,1} denotes the binary 
function which specifies the transmission of the first material grating. The latter serves 
to imprint a density modulation onto the beam, thus creating the required spatial 
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coherence downstream at the light grating. We characterize the material grating mask 
by the Fourier coefficients 

Aj := - / T 1 (x)e- 27rlJ:c/d dx, (25) 

« J-d/2 

where d is the common period of all gratings. 

In order to evaluate the interference pattern ([24"j) it is now convenient to define the 
coherent Talbot-Lau coefficients in terms of the Fourier coefficients ([T2"[) which describe 
the phase shift due to the optical dipole potential, 

Bm (£) := E h ^-rn eX P (" 2 ™ [j ~ j] t) ■ W 

jez 

The incoherent effect of the light grating is best accounted for through the characteristic 
coefficients associated to the Fourier coefficients of the periodic probability distribution 
([T5"[) . They are given by 

XAO ■■= E P t ] ex P (-27ri*0 , (27) 

kez 

and they serve to define the general Talbot-Lau coefficients, which include the effect of 
absorption, 



5 m (£) := E 5 «(0Xm-n 



(28) 

The factor ~ in f [251) reflects the fact a single photon has a momentum that equals one 
half of the grating momentum p d . As one expects, the convolution (1251) reduces to the 
coherent expression ( [261) if absorption can be neglected, i.e., for % m (£) = <5 mi0 . 

Inserting the Fourier expressions ([TBI and (|15|) into ([211) ^ ne integrations can now be 
carried out by retaining the resonant contributions. This yields the interference pattern 
in terms of the coefficients ([251) and (J25J, 

w 3 (x) = E A *Ae [Zj-) ex P ( 2 ™4) • ( 29 ) 

Here = d 2 /Adb denotes the Talbot length, which gives the characteristic length scale 
for near- field interference. 

One records the beam intensity behind the third grating as a function of the 
lateral position x s , S (x s ) oc J dxT 3 (x — x s ) w 3 (x). Since the first and third grating 
are identical in our experiments, T\ (x) = T 3 (x), the expected interference signal reads 

S (x s ) = £ (A*) 2 B 2e U±-\ exp (2ttz^) . (30) 



Classical description. Using the same general formalism as above the classical result 
is obtained by replacing the Wigner function by the classical phase space density 
and the kernel T co \ 1 {x,p) in ([241) by its classical counterpart ( [231) . The evaluation of 
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the corresponding moire-type density distribution suggests to introduce the classical 
coefficients 

C m (£) = I f d/2 dxexp (-2mm-\ exp ( -2m^-^-A . (31) 
dJ-d/2 v dJ \ Pd J 

They are the classical analogue of the Talbot-Lau coefficients (1261) . but clearly lacking an 
interference phase factor. Performing the same steps as above, the classical prediction 
for the signal behind the third grating thus assumes a form analogous to ( |30l) . 



Sd (x s ) = Yl ( A *rf ( £ f-) exp ( 2mi l[) ' (32) 

where Lt = d 2 mv z /h, with m the molecular mass and v z their longitudinal velocity. 
Note that Planck's constant appearing in Lt cancels against the one from pd = h/d 
showing up in (l3"Tj) ; it is kept here to maintain the close analogy with the quantum 
result. Like in (1281) . the possibility of photon absorption is accounted for in (1321) by a 
convolution with the characteristic coefficients (|27|) . 



(33) 



S.^. Closed expressions for the sinusoidal light grating 



The results obtained so far are valid for gratings with arbitrary eikonal phase shifts 
and momentum kick distributions. We now focus on the complex light grating of our 
experiment, as defined by the intensity distribution ([T]). Their special form will yield 
closed formulas for the Talbot-Lau coefficients ( 1281) and their classical analogues ( 1331 . 

The sinusoidal x-dependence of the phase shift flHJ) implies that the Fourier 
coefficients are determined by the integer Bessel functions, 

b m = {-iTe^jJ^Y (34) 



The summation for the coherent Talbot-Lau coefficient ( 1261) can be carried out by means 
of Graf's addition theorem for Bessel functions EHI- This leads to 



B m (0 = -Jm (-00 sin «)) , (35) 

indicating that all the Talbot-Lau coefficients are real. 

It is instructive to compare this to the corresponding coefficients ( I3T1) of the classical 
formulation. It follows immediately from ([H]) that they are given by 

C m (0 = Jra (-7T^) ■ (36) 

Comparing the the quantum expression ( |35l) and the classical one (1361) . we see that both 
coefficients assume the same limiting form if the Talbot parameter £ is much smaller 
than unity. They do however strongly deviate for £ ^ 1. The distinguishing quality 
of the quantum wave coefficients ( 1351) is their periodicity in £, which gives rise to the 
characteristic Talbot-Lau recurrences. Classical particles show no such recurrences since 
their coefficients ( |36l) exhibit no periodicity in £. 
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We move on to evaluate the characteristic coefficients ( 1271) . There is no obvious 
way to express the Fourier coefficients P% from Eq. (ITS]) in closed form. However, 
the coefficients Xm{0 can be expressed as the Fourier transform of the characteristic 
function (T2"U1) with respect to position, 

r i/2 

Xm (0 = / dr exp (-2mmr) X (£, rd) . (37) 

.7-1/2 

Due to the sinusoidal position dependence of the mean photon number (jHJ) the 
integration can be carried out, yielding a modified Bessel function, 

Xm (0 = exp (-no sin 2 (7r£)) ("o sin 2 (7r£)) . (38) 

The Talbot-Lau coefficients in the presence of absorption can now be obtained by 
performing the summation in (j28p . This can be done using an addition theorem for 
mixtures of regular and modified Bessel functions, which can be derived from Graf's 
addition theorem. It reads, for u, v e R, u ^ v , 

/ _ \ n / 2 

\v + uj k ^ 

and yields a real number also for \v\ < \u\ since J n (iu) = i n I n (u). 

Using (1351) and (1371) . and noting J_ n (z) = J n (—z), we thus obtain the general 
coefficients of the Kapitza-Dirac Talbot-Lau interferometer, which incorporate the effect 
of photon absorption. They are given by 

' Ccoh(0-Cabs(Q \ m/2 
Xcoh(0+Cab S (0/ 

x J m (-sgn [C abs (0 + Ccoh(01 - CL(o) • (40) 

Here, the coherent diffraction effect of the dipole force is described by the function 

Ccoh(O = 0osin(7rO (41) 
and the incoherent effect of absorption is accounted for by 

Cabs(0=^osin 2 (^) • (42) 

Classical description. The coefficients for the classical motion (l3"3"j) can be obtained the 
same way starting from ( |36l) and ( TSTl) . Given the relation between the classical and the 
quantum coefficients (IHUj) and (1531) . it is not surprising that the C m (£) assume a similar 
form as the B m (£) in (HOI) . The only difference is that Ccoh(0 is replaced by 

Cci(£) = 0o7r£, (43) 
which lacks the periodicity in the Talbot parameter £ shown by (|4ip . 



S m (£) = exp (-Cabs(O) 
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2. 5. Discussion of the theoretical results 

2.5.1. The fringe visibility Using the above results it is now easy to calculate the 
expected quantum interference pattern ( 130]) and the corresponding classical prediction 
( 13*2]) . However, for the parmeters of our experiment the patterns are well described by 
a sine curve so that it is sufficient to characterize the experimentaly observed pattern 
by the contrast of a sinusoidal fit. This sinusoidal fringe visibility can be calculated as 
the ratio of the first two Fourier coefficients of the fringe pattern, V = 2\S\/ Sq\. 

We denote by / the open fraction (i.e., the ratio between the single slit width and 
the grating period) of the first and the third grating, so that the grating coefficients 
( I25p are given by At = / sine (faf). The quantum fringe visibility then takes the form 



Similarly, the fringe pattern expected form classically moving molecules has a 



function of coherent diffraction Ccoh(£) by the classical function Cci(C) defined in Eq. (113]) . 

2.5.2. Comparing the classical and quantum predictions. Let us now see how the 
quantum interference pattern differs from the fringe pattern expected from classically 
evolving particles. Figure [2J compares the corresponding visibilities V qm and V c i as 
one varies the de Broglie wave length or velocity. The latter is specified by L/Lt = 
L/d 2 x AdB both in the quantum and the classical case. For the material gratings we 
assume an open fraction of / = 0.42 , like in our experiment. 

The most important feature of the quantum result (top row in Fig. [2j is that the 
visibilities are generally much greater than the classical calculation (bottom row). They 
also display more structure if the strength of the dipol potential increases (left column 
to right column), a consequence of the intricate near field interference process. At 
first sight, it may seem surprising that a fringe visibility would be observed at all in 
this setup if the molecules were moving as classical particles. This is due to a moire- 
type effect, where the light grating acts as a periodic structure of lenses focusing the 
classical trajctories. Note that these classical visbilities are systematically suppressed 
in the 'quantum regime' L/L^ ^> 1. They coincide with the quantum result only in the 
'classical limit' L/L^ — > of a vanishing de Broglie wave length. 

One also observes in Fig. [2J that the visibility peaks are affected rather differently 
by the possibility of photon absorption. Close to even multiples of L/L-r the classical 
contrast remains essentially unaffected by absorption, while the quantum visibility 
vanishes identically at all integer multiples of L/Lt in the absence of absorption. It 
is a curious result of our theory that a certain fringe pattern can be obsered even if the 
dipole force can be neglected compared to photon absorption, O ~> 0. The quantum 
and the classical predictions coincide in this case, and one expects a sinusoidal visibility 




(44) 



with the coefficient B 2 given by (14*0]) . 




obtained from ([44]) . if we replace the 
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Figure 2. Fringe visibilities as obtained from the quantum (top) and classical 
(bottom) description of the molecular motion in a KDTLI. The abscissa value L/Lt 
is proportional to the molecular de Broglie wave length; it scales as v~ x , also in 
the classical case. The material gratings are chosen to have an open fraction of 
/ = 0.42 13J. The maximal optical phase shift increases from 4>o = 3 (left panel), 
<f>o = 5 (middle) to cf>o — 7 (right). Each line in each panel describes the fringe visibility 
for a specific maximal mean number of absorbed photons no, which is given, from the 
top to the bottom lines, by no/(f>o = 0%, 10%, 25%, and 50%. The quantum and the 
classical fringe visibilities differ markedly both in their value and in their functional 
dependence, except for the limit of high velocities, i.e. small de Broglie wave lengths 
L/Lt — > 0. The visibility peaks are strongly affected by photon absorption. 



given by 



V a bs = 2 sine 2 (itf) exp 



-c 



abs 



L 



abs 



L 



(45) 



It is greatest if the maximum mean number of absorbed photons equals n = 4.65. The 
visibility then amounts to max(V a bs) = 0.236 x sine 2 (nf), i.e., a value of 24% cannot 
be exceeded by this effect. 

Finally, Fig. [3] shows for a fixed value of L/L^ how the visibilities depend on 
the molecular properties, which are summarized in the dipole force phase 0o an d the 
absorption number n . We choose L/L T = 8.5, which corresponds to C 60 fullerenes at 
a velocity of 97m/s. One observes that at fixed velocity the parameter dependence 
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is less complicated than the wave length dependence of Fig. [2j Since the molecular 
velocity is easy to control this implies that KDTL-Interference can be used to measure 
the molecular polarizability and the absorption cross section by varying the intensity of 
the light grating [EH 123 ED]- 

So far, the molecular beam was assumed to be characterized by a fixed longitudinal 
velocity v z . The case of a finite velocity spread is easily incorporated by averaging 
the interference patterns of the monochromatic theory with the measured velocity 
distribution in the beam. This also applies to the sinusoidal visibilities, since the zeroth 
Fourier component of the interference pattern is independent of the velocity. 




Number of absorbed photons n Standing light wave laser power (W) 



Figure 3. Interference visibility for a fixed value of L/Lt — 8.5 and / = 0.42 for Cgo, 
which corresponds to a velocity of 97 m/s. Part a) shows the quantum interference 
visibility as a function of $o an d no. Panels b) - d) display cuts through this figure, 
along with the smaller classical fringe visibility (solid lines), b) is the cut along no = 0.5, 
c) shows the cut for fixed <t>o — 2.7, and d) is a cut along the surface shown in a), which 
corresponds to a linear increase in the power of the light grating according to equations 
(jl]) and 0. Also shown in d) is the effect of a velocity distribution Av/v — 10% on 
the quantum (dashed dotted curve) and classical (dash double dotted curve) visibility. 

3. Verification of the model using fullerenes and fiuorofullerenes 

In a first experiment a mixture of the fullerenes Ceo (720 amu) and C70 (840 amu) was 
co-sublimated in a thermal source (see Fig. [T]) at a temperature of 910 K. By alternating 
the setting of the detecting quadrupole mass spectrometer between the masses of the 
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Laser Power (W) 

Figure 4. Power dependence of the interference contrast for the fullerenes Cgo 
(filled circles) and C70 (open circles). The points represent the weighted average 
of three consecutive measurements, the error bars depict one standard deviation. 
Bold lines display the quantum expectations, obtained by weighting Eq. ([44]) with the 
experimentally observed velocity distribution and fitting a and 00 • The corresponding 
classical expectations are shown as thin lines. The solid lines identify the theoretical 
predictions for Ceo while the dashed lines correspond to C70. 

two molecules, the interference patterns of either species were recorded, one after the 
other, before the laser power was shifted to the next point in the series. This way, 
we assured that the standing light wave conditions were the same for both molecules. 
The result is shown in Figure HI Each data point represents the weighted average of 
three consecutive measurements, where the confidence intervals of the individual sine 
fits serve as weights. The single patterns were recorded over four full sine periods with 
ten points per period and two seconds integration time per point. The mean count rates 
per second amounted to 740 for Ceo and 505 for C70. 

The velocity of the particles was selected by limiting their path to the associated 
free flight parabola in the Earth's gravitational field [40J. This yields an approximately 
Gaussian velocity distribution. The mean velocities in the experiment were determined 
to be 202 m/s for Ceo and 194 m/s for C70 with velocity spreads of 27% and 25%, 
respectively (Av/v, FWHM). The interferometer setup is characterized by a grating 
period of <i=266nm (corresponding to a laser wave length of 532 nm), a grating 
separation of L=105mm, a molecular beam height of about 200 /mi, and a beam width 
of about 1mm. The slit widths are assumed to be 85 nm in Gl and 110 nm in G3. 

Both visibility curves reproduce the quantum expectations accurately (bold lines, 
Eq. (130|) ). while being in distinctive disagreement with the classical prediction (thin lines, 
Eq. (132jl ). We emphasize that the present result signifies a noticeable improvement over 
previous measurements [IB]. We attribute the enhanced interference contrast mainly 
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to further improvement of the highly critical adjustment of the machine as outlined in 
Section 0J 

From the perfect accordance with the theory we also deduce that decoherence due 
to collisions with particles of the background gas is negligible in our current experiments. 
Following [JT] we estimate an effective cross section for collisions of <7 e fr = 4.2 x 10~ 17 m 2 
for C70. The experiment is conducted under pressures below 10~ 8 mbar, which results 
in a mean free path of more than 17 meters. This corresponds to a reduction of the 
effective visibility of V e s > 0.98 Vo- Since the effective collisional cross section is mainly 
governed by the polarizability of the molecule rather than its geometrical size, an even 
smaller reduction of the visibility can be expected for the other species discussed in this 
article. 

We also observe, and again in good agreement with our model, that the fringe 
contrast of the more absorptive C70 decays significantly more rapidly than the contrast 
of Ceo when we increase the laser power. The increasing number of absorbed photons 
fills in the interference minima with shifted interference curves, thus washing out the 
accumulated interference pattern. 

It may come as a surprise that, in spite of the higher absorption cross section, C70 
exhibits actually a higher interference contrast than Cqo at lower laser powers. This can 
be explained by the optical polarizability which, according to our present measurement, 
amounts to ctAc — 114A for C70 and is thus 31% higher than for Ceo- This results in a 
larger phase shift in the optical grating and thus leads to a shift of the entire curve to 
the left in Fig. @J 

The experiment was repeated in two additional and separate runs with the 
fluorofullerenes C60F36 and C60F48. The C60F36 sample was synthesized by the method 
described in [42] with a compositional purity of > 95% as determined by mass spectrum 
analysis. It contains two major isomers of C3 and Ci symmetry in the approximate 
ratio 2:1, as well as one minor isomer of T symmetry (ca. 5%) [43]. The three 
isomers of C60F36 show very similar polarizabilities, ranging between 62A and 65A, 
according to calculations performed with Gaussian |44j . Depending on the specific 
isomer, C 60 F 36 may posses a dipole moment of up 1.2 Debye. The sample of C 60 F 48 was 
purchased from Prof. L. Sidorov, Moscow. The synthesis and characterization was done 
according to [45]. C60F48 is formed predominantly as a single isomer of D3 symmetry, 
with the minor isomer of Sq symmetry comprising ca. 5% |43J. For both isomers of 
C60F48 the simulations yield virtually identical polarizabilities and no dipole moment. 
Both fluorofullerene samples were produced before the year 2003 but remained intact 
compounds over this period, as proven by mass spectra. 

It is noteworthy that an earlier experiment in a pure Talbot-Lau configuration 
succeeded already with C 60 F 48 but at limited contrast [SI] . The present setup however 
substantially outperforms its predecessor: The more sophisticated KDTLI scheme, an 
improved count rate, and a better vibration insulation with respect to the earlier 
experiment allowed for the first time to achieve the full expected quantum contrast 
for both C 60 F 36 and C 60 F 48 , as shown in Fig. El The larger error bars with respect to 
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Figure 5. Power dependence of the interference contrast for the fluorofullerencs 
CeoF36 (filled circles) and C60F48 (open circles). The circles represent the weighted 
average of three consecutive measurements, the error bars depict one standard 
deviation of the shot noise. Bold lines display the quantum expectations. The values 
used for the parameters a and 00 correspond to a best fit and are depicted in Tab. [3] 
Classical expectations are shown as thin lines. Solid lines correspond to C60F36, dashed 
lines to C60F48. 

the C60-C70 measurement are mainly a consequence of the lower count rates of only 60 
and 75 per second for C 60 F 36 and C 60 F 48 , respectively The temperature was kept a 
590 K for both molecules. 

The recorded velocities for C 6 oF36 and C 60 F 48 were 130 m/s and 116 m/s with 
velocity spreads of 16% and 18%, respectively. Although smaller velocities tend to make 
the experiment more susceptible to vibrations, drifts and misalignment, no significant 
drop of the measured visibility below the theoretical expectation was observed. 

In Tab. [3] we compare the optical properties of all four particles that were 
extracted from a best fit of the quantum curves to the experimental data. All values 
are in good agreement with the parameters determined in independent experiments 
[46, 29J and with molecular simulations carried out using Gaussian. The remarkably 
small statistical errors indicate that our method offers the capacity for high precision 
metrology experiments with heavy molecules. The accuracy is, however, currently 
limited by the systematic errors which are primarily governed by the accuracy of the 
measurement of the power (± 5%) and the waist (± 10%) of the diffracting laser beam. 

The decrease in polarizability from Cgo to C60F48 is in good accordance with the 
observation that fluorinated molecules in general show a reduced polarizability-to-mass 
ratio and correspondingly lower inter-molecular binding, lower sublimation enthalpies 
and higher vapor pressures at a given temperature [50]. It is also important to see that 
the fluorine shell reduces the absorption cross section at the wavelength of the diffracting 
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Molecule 


^60 


C70 


C60F36 


C60F48 


a ab J10- 22 m 2 l 


2.8 ± 0.3 ± 0.3 


24.9±1.1±2.7 


< 0.6 


< 0.5 


a opt [A 3 } (Exp.) 


87.1±0.5±9.7 


114.2 ± 0.9 ± 
12.7 


60.3±1.0±6.7 


60.1±0.8±6.7 


a sta t[A 3 l (Lit.) 


88.9 ± 0.9 ± 
5.1 gS] 


108.5 ± 2.0 ± 

6.2 [ig 


62-65 [34] 


63 [S] 


a [A 3 ] (Lit.) 


89.2 07J 


109.2 071 

L — —i 






a [A 3 ] (Lit.) 


90 [IE] 


118.4 08] 






a [A 3 ] (Lit.) 


98.2 [49J 


122.6 09] 







Table 1. Molecular parameters as derived from the best fit of the theory including 
statistical and systematic errors. The data are in very good agreement with the 
values provided by literature, where published, or with simulations performed with 
Gaussian [44]. The simulation yields slightly different values values for the three 
conformers of C6oF36, ranging between 62 A 3 and 65 A 3 . Note that the published 
and calculated values represent static polarizabilities while our experiment yields the 
optical polarizability at the laser wavelength of 532 nm. 

laser beam to a negligible value. Our measurement thus allows us to extract information 
about the effect of fluorination on the electronic properties of fullerenes. 

4. Alignment requirements for precision experiments 

Matter wave interferometry with large molecules operates with de Broglie wavelengths in 
the range of a few picometers and grating periods as small as a few hundred nanometers. 
As a result of that, the interferometer alignment has to be considered carefully. The 
following section is therefore devoted to a short assessment of the constraints on the 
experimental precision. 

4-1. Equality of grating periods 

If the first two gratings have only slightly different lattice periods, the interference 
pattern spacing will not match the period of the third mask and the contrast will be 
reduced. A period mismatch as small as one per mille leads already to half a fringe shift 
between slit one and five hundred. 

In practice, all grating periods must be, on average, equal to better than 0.05 nm, 
i.e about the diameter of a hydrogen atom. This condition enters both the choice of 
the grating manufacturing process and the alignment of the yaw angle for all gratings. 
Modern photo-lithography and etching procedures allow to reach this level of precision. 
The gratings for our experiments were produced by Dr. Tim Savas at MIT and 'nm 2 ' Inc, 
Cambdridge, Massachusetts, and independently checked by Ibsen photonics, Denmark. 
The gratings were fabricated to be 0.3 nm wider than the period of the standing light 
wave in order to allow for later yaw adjustments. 
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4-2. Transverse grating shifts and grating roll angles 

The lateral position Axi of all gratings relative to each other determines the final location 
of the fringe pattern. For Talbot-Lau interferometry we can define a phase of the 
near-sinusoidal interferogram which is determined by the relative shift of the molecular 
density maxima with regard to the openings of the third grating. For geometrical reasons 
this phase is determined by (e.g. [5"T] ) 



with kd = 2n/d. In a symmetrical setup such as ours, with L\ = L 2 , the prefactor 
of the first and third grating must be equal while the second grating's shift enters 
twice. Molecules passing the gratings at different heights y\ will effectively see different 
transverse grating shifts Ax; = aiy\ if element i is rolled by the angle a.\. Although a 
simple phase shift will not reduce the fringe visibility, an integration over many height- 
dependent shifts does. If we neglect gravity, we can assume the second grating to define 
the reference angle, a 2 = 0, and study the influence of rolling Gi by ct\ and G3 around 
a 3 . In the paraxial limit the total signal Six) behind the third grating is then 



where the transmission function S is parametrized as a sine wave with visibility V and 
height dependent phase shifts x\ and x 3 . We denote a[ = aiL 2 (Li + L 2 + L 3 ) / (LiL tot ) + 
a 3 L 3 /L tot , a 3 = ol x LqL 2 I (L 1 L tot ) + a 3 (L + L x + L 2 ) / L tot with L tot = (L + L 1 + L 2 + L 3 ) 
and all distances as shown in Figure El 

From Eq. ( T4*8l) we see that rolling Gi and G3 reduces the fringe visibility in a 
sinc-shaped functional dependence. The alignment becomes increasingly important for 
smaller grating constants and more extended molecular beams. In our experiment, the 
first zero of the sinc-curve appears for a roll angle of around 0.65 mrad. 

In addition to shifting the phase, rolling also affects the effective grating constant. 
If one of the gratings is rolled with respect to the others its projected period increases 
by d cf[ = dj cos(a) ~ d + da 2 /2 + 0(a 4 ) (see Fig. E])- If we require the relative period 
change not to exceed (d e g — d)/d= 10~ 4 , the roll angle has to be kept aligned to within 
10 mrad. 

4-3. Longitudinal grating shifts 

The semiclassical picture in Fig. Eh shows that the interference pattern is also blurred 
when the third mask is moved relative to the second grating by the distance ±AL. For 
symmetry reasons, the same is true for a movement of GI. The contrast is severely 
reduced when the blur D = d/2 is as wide as half a grating period. We see that this 



4> = k d (Axi - 2Ax 2 + Ax 3 ) 



(46) 



h s /2 h d /2 




(47) 



-h 3 /2 -h d /2 

So(l + Vsin(x) sine (kdhga^) sine (kdhda' 3 )). 



(48) 
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Figure 6. Illustration of the alignment considerations: a) An imbalance of grating 
separations leads to blurred interferences; b) and c) Illustration of the grating's motion 
during roll, pitch and yaw. d) Distances relevant for the assessment of grating roll: 
molecules can pass the gratings at different altitudes and different transverse positions 
if the gratings are rolled with regard to each other. The height of the first velocity 
selection slit is h$ = 150 /im while the third one is restricted to = 200 /xm. The 
distance between the source and G\ is Lq — 150 cm while the spacing between G3 and 
the detector is L3 = 250 mm. The grating separation amounts to L\ = L2 = 105 mm. 
The beam width of about 1 mm is sufficient to illuminate nearly 4000 grating openings. 
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condition is met when 

D d/2 Nd 



(49) 



AL AL 2L ' 

where N is the number of grating slits illuminated by the molecular beam. We thus 
derive the length-balance criterion 
AL 1 

- < w < 50 > 

When 4000 lines are illuminated, as in our experiments, we have to balance the distances 
to better than 25 fim. This is already comparable to the waist of the diffracting laser 
beam. This intuitive condition is consistent with a complementary and more rigorous 
treatment using Wigner functions [12] . 



4-4- Grating pitch 

The effect of forward or backward tilting a single grating, i.e. to add a pitch, is to 
introduce a height-dependent imbalance in the grating separation. The pitch must 
be compatible with the requirement of Eq. (I50p . If the beam height is h and the 
forward pitch is measured by the angle 9 then the arm lengths are balanced as long as 
Az = h9 <C L/N. For our experiment with h = 100 //m, L = 105 mm and N = 4000, 
this corresponds to 250 mrad. This condition can easily be met. 



4-5. Grating yaw 

To first order, the argument for grating pitch also holds for grating yaw. However, yaw 
also changes the slit's effective period as well as their effective open width. Under a yaw 
of angle <ft the grating constant shrinks like d y = cos(<f))d « d — (p 2 d/2. Similarly to the 
roll-related period change we derive the condition < 10 mrad. 

Finally, we have to include the effective reduction of the open slit width if the grating 
is turned: The openings shrink because of the finite wall thickness b to a e g = a — b tan (p. 
A reduction of the open fraction in the first or third grating is important as it tends to 
increase the fringe visibility while decreasing the count rate at the same time. A variation 
by 10% in visibility is already rather clearly noticeable, leading to the constraint: 
0.1 > (a e ff — a) I a = — (6/a)tan0. For a = 90 nm and b = 190 nm we thus find a 
limit of = 47 mrad for the maximally allowed yaw angle in our experiments. 

For the second grating a different reasoning applies: The optical grating is about 
20 yum thick along the direction of the molecular beam, therefore we have to make 
sure that no molecule crosses a sizeable fraction of the standing light wave period 
transversely. Since the period is as small as 266 nm the condition imposes a limit 
on both the angle of incidence and the divergence angle of the molecular beam: 
(f) < 0.1 x 266/20, 000 ~ 1 mrad. 
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5. Conclusions 

Our experiment combines the virtues of near field interferometry with the advantages of 
optical manipulation. Compared to far-field diffraction of collimated beams, near-field 
interferometry provides higher signal throughput at the expense of increased alignment 
requirements. As discussed in this paper, the experimental challenges are non-negligible 
but manageable with reasonable effort. Optical phase gratings allow us to operate with 
highly polarizable molecules, which would otherwise acquire prohibitively large van der 
Waals phases by interacting with the walls of material gratings. 

The Kapitza-Dirac- Talbot-Lau interferometer is very well described using a Wigner 
function approach, which facilitates in particular the inclusion of grating transformations 
and momentum exchange in the interferometer. The model is in excellent accordance 
with the experiments, even for particles as complex as C60F36 and C60F48. Actually the 
KDTLI allowed us to observe much improved contrast compared to earlier pure Talbot 
Lau experiments with the same particles |31j. 

This instrument and its theoretical description are now a good basis to continue 
the quest for the ultimate mass and complexity limits of matter-wave interferometry. 
In addition, it has proven to be useful for a number of relevant measurements of 
molecular properties, such as optical polarizabilities and absolute optical absorption 
cross sections [291 [30J, [52] • Interestingly the rich internal structure of complex molecules, 
including electric dipole moments or magnetic moments, structural properties etc. can 
be ignored unless they either modify the optical polarizability or unless we introduce 
additional field gradients or collisions jH] which allow to couple to these properties 
separately [35] . 

Acknowledgment 

We thank the Austrian science foundation FWF for support within the projects Z1^9 
Wittgenstein and the doctoral program W1210 CoQuS. We also acknowledge financial 
support by the MIME project within the ESF Eurocore EuroQUASAR program. K. H. 
is supported by the DFG within the Emmy Noether program. 

[1] Tonomura A 1987 Rev. Mod. Phys. 59 639-669 

[2] Batclaan H 2007 Rev. Mod. Phys. 79(3) 929-941 

[3] Sonncntag P and Hasselbach F 2007 Phys. Rev. Lett. 98 200402 

[4] Rauch H and Werner A 2000 Neutron Interferometry: Lessons in Experimental Quantum 

Mechanics (Oxford Univ. Press) 
[5] Berman P R 1997 Atom Interferometry (New York: Acad. Press) 
[6] Cronin A D, Schmiedmayer J and Pritchard D E 2008 Rev. Mod. Phys. 

[7] Arndt M, Nairz O, Voss-Andreae J, Keller C, der Zouw G V and Zeilinger A 1999 Nature 401 
680-682 

[8] Arndt M and Hornberger K 2009 in Proceedings of the International School of Physics on 
"Quantum Coherence in Solid State Physics" , Course CLXXI, IOP, in press 

[9] Bruchl R, Guardiola R, Kalinin A, Kornilov O, Navarro J, Savas T and Toennies J P 2004 Phys. 
Rev. Lett 92 185301 



Theory and experimental verification of Kapitza-Dirac- Talbot-Lau interferometry 24 



[10] Brezger B, Hackermiiller L, Uttcnthalcr S, Pctschinka J, Arndt M and Zeilinger A 2002 Phys. Rev. 
Lett. 88 100404 

[11] Griscnti R E, Schollkopf W, Toennies J P, Hegerfeldt G C and Kohlcr T 1999 Phys. Rev. Lett. 83 
1755 

[12] Nimmrichtcr S and Hornbcrger K 2008 Phys. Rev. A 78 023612 

[13] Gerlich S, Hackermiiller L, Hornberger K, Stibor A, Ulbricht H, Gring M, Goldfarb F, Savas T, 

Miiri M, Mayor M and Arndt M 2007 Nature Physics 3 711 - 715 
[14] Kapitza P L and Dirac P A M 1933 Proc. Camb. Philos. Soc. 29 297 - 300 
[15] Frcimund D L, Aflatooni K and Batelaan H 2001 Nature 413 142 - 143 

[16] Moskowitz P E, Gould P L, Atlas S R and Pritchard D E 1983 Phys. Rev. Lett. 51 370 - 373 
[17] Gould P L, Ruff G A and Pritchard D E 1986 Phys. Rev. Lett. 56 827 - 830 
[18] Martin P J, Oldaker B G, Miklich A H and Pritchard D E 1988 Phys. Rev. Lett. 60 515-518 
[19] Steane A, Szriftgiser P, Desbiolles P and Dalibard J 1995 Phys. Rev. Lett. 74 4972 - 4975 
[20] Rasel E M, Oberthaler M K, Batelaan H, Schmiedmayer J and Zeilinger A 1995 Phys. Rev. Lett. 
75 2633 - 2637 

[21] Giltncr D M, McGowan R W and Lee S A 1995 Phys. Rev. Lett. 75 2638 - 2641 
[22] Cahn S B, Kumarakrishnan A, Shim U, Sleator T, Berman P R and Dubetsky B 1997 Phys. Rev. 
Lett. 79 784-787 

[23] Turlapov A, Tonyushkin A, Sleator T, 2005 Phys. Rev. A 71 043612 
[24] Tonyushkin A, Sleator T, 2006 Phys. Rev. A 74 053615 

[25] Ovchinnikov Y, Mullcr J, Doery M, Vredenbregt E, Helmerson K, Rolston S and Phillips W 1999 

Phys. Rev. Lett. 83 284- 287 
[26] Deng L, Haglcy E W, Wen J, Trippenbach M, Band Y, Julienne P S, Simsarian J E, Helmerson 

K, Rolston S L and Phillips W D 1999 Nature 398 218 
[27] Nairz O, Brezger B, Arndt M and Zeilinger A 2001 Phys. Rev. Lett. 87 160401 
[28] Brezger B, Arndt M and Zeilinger A 2003 J. Opt. B: Quantum Semiclass. Opt. 5 S82-S89 
[29] Hackermiiller L, Hornberger K, Gerlich S, Gring M, Ulbricht H and Arndt M 2007 Applied Physics 

B 89 469-473 

[30] Gerlich S, Gring M, Ulbricht H, Hornberger K, Tiixen J, Mayor M and Arndt M 2008 Angew. 

Chem. Int. Ed. 47 6195 - 6198 
[31] Hackermiiller L, Uttcnthalcr S, Hornberger K, Reiger E, Brezger B, Zeilinger A and Arndt M 2003 

Phys. Rev. Lett. 91 90408 
[32] Patorski K 1989 in E Wolf, ed, Progress in Optics XXVII (Amsterdam: Elsevier Science Publishers 

B. V.) pp 2-108 
[33] Clauser J F and Li S 1994 Phys. Rev. A 49 R2213 

[34] Pfau T, Spalter S, Kurtsiefcr C, Ekstrom C and Mlynck J 1994 Phys. Rev. Lett. 73 1223 1226 
[35] Chapman M S, Hammond T D, Lenef A, Schmiedmayer J, Rubenstein R A, Smith E and Pritchard 

D E 1995 Phys. Rev. Lett. 75 3783-3787 
[36] Hornberger K, Sipe J E and Arndt M 2004 Phys. Rev. A 70(5) 53608 
[37] Hillcry M, O'Conncll R F, Scully M O and Wigner E P 1984 Phys. Rep. 106 121 
[38] Englcrt B, Miniatura C and Baudon J 1994 J. Phys. II France 4 2043-2059 
[39] Erdelyi et al. Tables of Integral Transforms vol I (McGraw-Hill, New York 1954) 
[40] Arndt M, Nairz O, Petschinka J and Zeilinger A 2001 C.R. Acad. Sci. Paris t. 2, Serie IV, 1-5 
[41] Hornberger K, Uttenthaler S, Brezger B, Hackermiiller L, Arndt M and Zeilinger A 2003 Phys. 

Rev. Lett. 90 160401 

[42] Boltalina O V, Borschevskii A Y, Sidorov L N, Street J M, Taylor R 1996 Chem. Commun. 4 
529-530 

[43] Popov A A, Senyavin V M, Boltalina O V, Scppclt K, Spandl J, Feigerle C S, Compton R N 2006 

J. Phys. Chem. A 110(28) 8645-8652 
[44] Gaussian 03W, V6.0, Frisch M et al., Gaussian, Inc., Pittsburgh PA, 2003 
[45] Boltalina O V and Galeva N A 2000 Uspekhi Khimii 69(7) 661-674 



Theory and experimental verification of Kapitza-Dirac- Talbot-Lau interferometry 25 



[46] Berninger M, Stefanov A, Deachapunya S and Arndt M 2007 Phys. Rev. A 76 013607 
[47] Guha S, Menendez J,Page J, and Adams G 1996 Phys. Rev. B 53 13106 

[48] Eklund P, Rao A, Wang Y, Zhou P, Wang K, Holden J, Dresselhaus M, Dresselhaus G 1995 Thin 

Solid Films 257 211 
[49] Sohmcn E, Fink J, Kratschmer W 1992 Z. Phys. B 86 87 

[50] Boltalina O V, Markov V Y, Borschcvskii A Y, Galeva N A, Sidorov L N, Gigli G, Balducci G 

1999 J. Phys. Chem. B 103(19) 3828-3832 
[51] Oberthaler M K, Bernet S, Rasel E M, Schmicdmayer J and Zeilinger A 1996 Phys. Rev. A 54 

3165-3176 

[52] Nimmrichter S, Hornberger H, Ulbricht H, and Arndt M 2008 Phys. Rev. A 78 063607 



